Time evolution and matter wave interference in Fermi condensates 
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We discuss matter wave interference of Fermi condensates in strongly coupled s-wave and p-wave 
channels and the time evolution of a single cloud upon release from trap. In s-wave systems, where 
the order parameter is a complex scalar, we find that the interference patterns depend on the relative 
phase of the order parameters of the condensates. In p-wave systems involving the mixture of two- 
hyperfine states, we show that the interference pattern exhibits a polarization effect depending on 
the relative orientation of the two vector order parameters. However, in p-wave systems involving 
a single hyperfine state, we show that this angular effect reduces to an overall phase difference 
between the two interfering clouds, similar to s-wave. Lastly, we also point out that p-wave Fermi 
condensates exhibit an anisotropic expansion, reflecting the spatial anisotropy of the underlying 
interaction between fermions and the orbital nature of the vector order parameter. 

PACS numbers: 03.75.Ss, 03.75.-b, 05.30.Fk 



I 

o 
o 



^ : 

in 
cn 
o 

o 



I 

o 
o 



X 



Matter-wave interference is a very powerful tool to 
study quantum phase coherence between atomic Bose 
Einstein condensates (BEC) P, 0, and spatial quan- 
tum noise of bosons in optical lattices Similar 
techniques can also be applied to study Fermi conden- 
sates, ^, '3, 0, H, where superfluidity can be tuned 
from the BCS to the BEC regime. These experiments 
may reveal that the time dynamics in the BCS regime 
is overdamped (large Cooper pairs can decay into two 
atoms), while in the BEC regime it is essentially un- 
damped (tightly bound molecules are stable) [ifl 
Matter-wave interference experiments of s-wave Fermi 
condensates may be readily performed, since stable con- 
densates already exist. For s-wave Fermi condensates 
in the BEC regime quantum interference effects are ex- 
pected to be similar to those of atomic Bose condensates, 
and the interference pattern should depend essentially on 
the phase difference of the order parameters between two 
interfering clouds. 

In contrast, it is more interesting to study interfer- 
ence effects in p-wave superfluids because of the vector 
nature of the order parameter. Many groups have re- 
ported some progress towards the formation of p-wave 
Fermi condensates in single clouds 0, 0, and in 
optical lattices where p-wave Feshbach resonances 
have been observed. For Feshbach resonances currently 
tried in single clouds atom losses have been significant, 
and the realization of stable p-wave condensates has not 
been achieved yet. However, other unexplored Feshbach 
resonances in single clouds may show less dramatic two- 
body dipolar or three-body losses as observed in optical 
lattices 15]. Thus, these experimental difficulties may be 
surpassed in the immediate future. 

We discuss in this manuscript time dynamics and 
matter-wave interference of s-wave and p-wave Fermi 
condensates. Our main results are as follows. While 
in atomic BEC and s-wave Fermi superfluids quantum 
interference patterns depend essentially on the relative 



phase of the two clouds, we find that in the p-wave 
Fermi superfluids there can also be a strong dependence 
on the relative angle between the two vector order pa- 
rameters, thus producing a polarization effect. Further- 
more, we show that p-wave Fermi condensates exhibit an 
anisotropic expansion, reflecting the spatial anisotropy 
of the underlying interaction between fermions and the 
orbital nature of the vector order parameter. 

We consider a system of fermions with mass m in two 
hyperfine states (pseudospins), labeled by greek indices 
a = 1,2. The Hamiltonian density is (with ft = fc^ = 1) 



2m 



+ C/cxt(r,t) 



where repeated greek indices indicates summation, "i/^J, 
(tAq) are creation (annihilation) operators of fermions 
in state a, and H{t) = J drH{r,t) is the Hamiltonian. 
The trapping potential is C/oxt(r, t) = J2j=x,y.z ^oi^Vj/'^^ 
where Ljj{t<Q)= LUj,o are constants. 

The generating functional for non-equilibrium pro- 
cesses associated with H{t) is 

Z{t) = TrU\t,to)exp [-(3{H{to) - iic.N^)]U{t,to), 

where is the number operator for fermions of type 
a, /ic( is the corresponding chemical potential, and 
(3 = l/T is the inverse temperature. Here, li{t,to) = 
exp[— i H(t')dt'] is the time evolution operator. This 
expression implicitly implies that the system is in ther- 
mal equilibrium at any time to < Oi since Z{tQ) = 
Trexp[—/3{H{to) — ^aNa)]. By introducing a complex 
time T, the generating functional can be written as 

Z{t)= [ n^D[^Pl,^^]e-'^^'^'''^^-'''^'^+^^^'''h'''^^\ (1) 
Jbc 

where the boundary condition (BC) of the functional in- 
tegral is ipa{Y,to — i/T) = — ■0Q(r, ^o)- The quadratic 
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FIG. 1: Integration contour C used in Eq. (|5J. 



term is 

320^1^,) = jj^ j rfr^l(r,T)£„(r,r)V^„(r,r), (2) 

where the integration contour C is shown in Fig. 1. The 
one-particle operator Ca is defined as 



2m 



In what foUows, we discuss the p-wave case in detail 
and quote the more standard results for s-wave. To cal- 
culate the fourth order term 6*4, we write Vap-ys(^) — 
V{'x)Taj3^S) where in a triplet channel Ta^^s — ^ai3 ■ 
(v'l')^^. The pseudospin matrices {vj)afi = {icrjO'y)a0, 
where aj are Pauli matrices. This form of interac- 
tion implies that the interaction strength is the same 
for all pseudospin channels. (Deviations from this 
assumption will be discussed later). By introducing 
B'^(r, r',r) = VlCr, T)va/3'!/'^(r', r) and the correspond- 
ing auxiliary field dt(r,r',T), the interaction term S4 
can be decoupled via a Hubbard-Stratonovich trans- 
formation. We integrate out the fermions and write 
d(r, r',r) = D„(R, T)r/„(x), where 77„(x) are eigen- 
functions of the reduced two-body Hamiltonian 0.2 — 
— Vx/to + V{x.). Here, R = (r -I- r')/2 and x = r — r' are 
center-of-mass and relative coordinates, respectively. In 
strong coupling limit where the interaction is the largest 
energy scale in the problem (with typical strength much 
larger than temperature and trapping potential), the 
ground state of H2 has the same symmetry of V. For 
definiteness, we consider a pure p-wave interaction where 
the ground state is three-fold degenerate and labeled by 
v — x,y,z, with corresponding eigenfunctions rio(x)x^. 
Thus, at low temperatures where the higher energy states 
are not excited, we have d = D^(R, T)rjQ{x)x^. 

In this strong coupling limit. Cooper pairs are tightly 
bound molecules and their relative degrees of freedom 
can be integrated out leading to 



(3) 
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dr j dR|Dt(R,r)- ii:D(R,r 
[2|D(R,r)r-|D2(R,r)p]}, 



where DjiTl) = Dj^^(R.) is the order parameter, and 
the operator K = idr — 2f/ext(R, t) + V^/(4m) corre- 
sponds to the action of an ideal non-equilibrium gas of 



Bose particles with mass M = 2m. This action leads to 
equations of motion 



idtDj 



2M 



+ 2C/e,t(R, t) + 2go\T>f 
.90 CD--D)d]. 



(4) 



Notice that this expression is different from the time- 
dependent Gross-Pitaevsku (TDGP) equation for a vec- 
tor boson field. The difference comes from the last term, 
which describes a non-unitary complex order parameter 
of the underlying paired fermions. In contrast, the stan- 
dard TDGP equation for scalar bosons is obtained in the 
s-wave case 16]. 

Equation Q can be simplified to the TDGP form 
in two special cases. First, if two hyperfine states are 
equally populated and D is unitary, then it can be rep- 
resented by a real vector with an overall phase factor, 
leading to the equation of motion 



^ + 2f/e,t(R,t) + go|D| 



(5) 



Second, if only one hyperfine state is populated, then D 
is non-unitary and D = A(l, —i, 0), where A is a complex 
constant. Thus, the last term D • D in Eq. Q vanishes, 
and the equation of motion is identical to Eq. (|SJ, with 
go ~f 2go. In the following discussion, we confine our- 
selves to these two special cases. 

For definiteness, we assume that the trapped Fermi 
superfluid is released at time t = 0, i.e., J7oxt(i > 0) = 0. 
Thus, for t < 0, the system is described by 



AioDj(R) 
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2t/c,t(R)+.9|D| 



(6) 



where /io is the effective boson chemical potential, g — 
go for the unitary case corresponding to two hyperfine 
states, and g — 2gQ for a non-unitary case corresponding 
to a single hyperfine state. When the composite Boson 
interactions are dominant the Thomas-Fermi (TF) ap- 
proximation leads to |D(R,0)| = g-i/2[^o-2J7oxt(R)]^/^ 
for ^0 > 2C/e:Et(R), and |D(R,0)| = 0, otherwise. When 
this approximation fails the initial condition for the time 
evolution can be obtained by solving Eq. ^ numerically. 

For t > 0, we use the transformation Rj{t) = 
&_i(t)i?j(0), where bj{t) are scaling factors satisfying, [itI 



oj: 



A{t)b,{t) 



(7) 



with initial conditions bj{t < 0) = 1 for all j, and A{t) = 
bx{t)by{t)bz{t). The D vector can be written as 



i?,(R(0,t) = 



VW) 



,(R(0),i)exp [zS{R{t),t)], 
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which upon substitution into Eq. Q and Eq. © at 

t = leads to 



S{R{t),t) = So{t)+MY, 



hk{t) dt 



(8) 



For a cigar-shaped trapping potential with axial symme- 
try (e = Wz(0)/a;_L(0) <C 1 and uj^_=ujx ^ ujy), Dj{Ii,t) 
becomes 

D,(R,t)«^^^H[^^I?,(R,0), (9) 

where Rk = Rk/hk{t) are scaled coordinates, and A = 
w_l(0)< is the dimensionless time. The result for s-wave is 
formaly identical to that of Eq. © with the substitution 
Dj — > 5*, where 4" represents the scalar order parameter. 

The matter-wave interference of two spatially sepa- 
rated condensates (without tunneling) is described by 



$tot(R, t) = $l(R, t) + $r(R, t), 



(10) 



where $p cx i'^j Dj^pajay denotes the pair wavefunc- 
tion of Fermi condensates in the left [P = L(+)] or right 
[R(— )] traps. Here, a coordinate system is chosen such 
that the trap centers lie at (-VF/2,0,0) and (VF/2,0,0), 
where W is the distance between the traps. We consider 
the case of two identical cigar shaped Fermi condensates, 
with the weakest confinement along the z axis, just like 
the experiment in Bose systems 0. In this case, the time 
evolution of the two independent Fermi condensates is 
described by 

Z,,..(R..) = -'''"'™^-"l l>,(R±Wx72,0,. 

Thus, for any single run of the experiment, the particle 
density n(R, i) = |$tot(R, i)P is 



n(R,t) cx |DL(R,i)P + |DR(R,t)|^ 



+2Re 



T>^ (R + Wi/2, 0) ■ D(R - W±/2, 0) 



AX 



,(11) 



where the phase x(R, t) = S{Ii + VKx/2, t) - 5(R - 
Wx/2, t) -I- Xo, and xo is the initial relative phase of the 
two clouds. The result for s-wave is formaly identical to 
that of Eq. lfTT|l with the substitution Dj ^ 

In the case where both Fermi condensates are in uni- 
tary states, D is essentially a real vector with an overall 
phase, and n(R, t) shows an angular dependence con- 
trolled by the dot product in Eq. (fTT|l . When the two 
order parameters are parallel, this term is maximal and 
the interference pattern is most visible (Fig. 2). However 
if the D vectors are perpendicular, fringes are absent at 
all times (Fig. 3). Therefore, in the unitary case the 
existence and intensity of interference fringes are very 
sensitive to the relative orientation of the vector order 
parameters. 
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FIG. 2: Interference pattern versus dimensionless time A for 
unitary p-wave Fermi condensates in the BEG limit involving 
two hyperfine states (cross section view). Each condensate is 
initially trapped in a cigar-shaped potential with Dj^ || Dr. 
The patterns are similar for atomic scalar bosons, and for 
s-wave and single hyperfine state p-wave fermions. 
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FIG. 3: Same plots as in Fig. 2, but with D[ ± Dr. 



When only one hyperfine state is occupied, the D vec- 
tor has a fixed form A(l,z,0). Without this pseudospin 
space rotational degree of freedom, the angular depen- 
dence of the interference pattern disappears and fringes 
are present in all experimental realizations. This result 
is similar to the s-wave case where the order parameter 
is a complex scalar. 

Before concluding, we would like to mention two im- 



4 



portant points. First, away from the BEC limit the ki- 
netic energy term in Eq. Q is 



E 



thus acquiring a mass anisotropy M^^^ = M/cl^^, which 
is directly related to the anisotropy of the Ginzburg- 
Landau coherence length ■ Here, we use the ba- 
sis of spherical harmonics Yi^mg and define Dj^rne by 
dj(v,v') = ^i,™«(R-)'7o,m«(a;)yL,mf(x). For a weak 
trapping potential, the coefficient c]'^ becomes 



Si 



x(k) /3r(k) 



4£;2(k) 16£:(k) 



/32k^x(k)y(k) ^1 



^2 J 



(12) 



where 
X(k) 



£;(k) = ^i.k -I- C2,k, Ca,k = kVSrn - /^Q, 
= tanh(/3a,k/2) -f tanh(/36,k/2), and y(k) = 
sech^(/3^i_k/2) + sech^(/3^2,k/2). The symmetry function 
is defined by T/(k, k') '= / dxV^(x) exp[i(k-k') -x] = 
F0(fc)0(fc')Yi^m^(k)Yj^*^/ (k'), and the angular average 



^, = 1 dkfc,%yi,™,(k)yi:„,(k). 



When the order parameter is characterized by one Yi^mi 
with mi fixed, k is diagonal in both mi and i, hence 



Cj,miSij. This mass anisotropy reflects the higher 



angular momentum (p-wave) nature of the order param- 
eter for paired fermions, and it is completely absent in 
s-wave Fermi and atomic Bose condensates. 

This effective mass anisotropy has a non-trivial influ- 
ence on the time evolution of condensates after release 
from the trap. Making the scale transformation R'j — 



i I \/^j,int, we conclude that the 



cloud expansion predominantly occurs in the strongest 
confined direction in the scaled space. Thus, the time 
evolution of a Fermi condensate initially trapped in an 
axially symmetric potential with sizes = Ly < 
leads to an anisotropic cross section of the cloud at any 
time after release when ^ My (^a, ^ £,y)- In Fig. 4, 
we show the cloud anisotropy ratio as a function of the 
effective mass anisotropy ratio. 

Second, in the discussion above we assume a symmet- 
ric interaction in pseudospin space, i.e.. Vim, ^^1212 and 
V2222 are identical. However, experimental results for 
p-wave Feshbach resonances show a a small but finite 
separation in different channels (1J|]. Different interac- 
tions in pseudospin space do not affect dramatically our 
derivation or results. When the different interactions are 
absorbed into an effective D vector, an equivalent proce- 
dure leads to equations similar to Eqs.(5) and (6). 

In conclusion, we considered a Fermi condensate con- 
sisting of two hyperfine states with s-wave and p-wave 
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FIG. 4: Cloud anisotropy ratio tl = Lx/ Ly as a function of 
effective mass anisotropy ratio tm = M^/My at time A (Solid 
lines). Dashed line indicate the saturated behavior at A — > 00. 



interactions, and derived equation of motion for p-wave 
case in a vector boson representation near the BEC limit. 
We found that the quantum interference of two p-wave 
Fermi condensates has a polarization effect due to the 
vector nature of the order parameter. This effect is ab- 
sent in strongly coupled (BEC) s-wave Fermi superflu- 
ids, as well as in scalar atomic Bose systems. Further- 
more, we observed that different orbital symmetries of 
the vector order parameter produce anisotropic effective 
masses (coherence lengths). For a cigar-shaped cloud 
with an isotropic cross section, the cloud expansion be- 
comes anisotropic at any time after release from the trap 
when the effective masses are anisotropic. Thus, the or- 
bital symmetry of the order parameter for p-wave con- 
densates can be probed via single cloud expansions. 
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